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ABSTRACT 

We derive an analytical theory of the prestellar core initial mass function based on an exten- 
sion of the Press-Schechter statistical formalism applied in cosmology. Our approach relies on 
the general concept of the gravo-thermal and gravo-turbulent collapse of a molecular cloud, with 
a selection criterion based on the thermal or turbulent Jeans mass, which yields the derivation of 
the mass spectrum of self-gravitating objects in a quiescent or a turbulent environment. With the 
same formalism, simply by using a constant density threshold for the selection criterion, we also 
obtain the mass spectrum for the non self-gravitating clumps produced in supersonic flows. The 
mass spectrum of the self-gravitating cores reproduces very well the observed initial mass func- 
tion, from the high mass domain to the brown dwarf regime, and identifies the different underlying 
mechanisms responsible for its behaviour. The theory predicts that the shape of the IMF results 
from two competing contributions, namely a power-law at large scales and an exponential cut-off 
(lognormal form) centered around the characteristic mass for gravitational collapse condition. 
The cut-off is not specifically due to turbulence and exists already in the case of pure thermal 
collapse, provided that the underlying density field has a lognormal distribution. Whereas pure 
thermal collapse produces a power-law tail steeper than the Salpeter value, dN/dlogM oc 
with X ~ 1.35, this latter is recovered exactly for the (3D) value of the spectral index of the ve- 
locity power spectrum, n ~ 3.8, found in observations and in numerical simulations of isothermal 
supersonic turbulence. Indeed, the theory predicts that x — (n+ l)/(2n — 4) for self-gravitating 
structures and x = 2 — n'/3 for non self-gravitating structures, where n' is the power spectrum 
index of log p. We show that, whereas supersonic turbulence promotes the formation of both 
massive stars and brown dwarfs, it has an overall negative impact on star formation, decreasing 
the star formation efficiency. This theory provides a novel theoretical foundation to understand 
the origin of the IMF and to infer its behaviour in different environment, characterized by the 
local properties of the gas, from today Galactic conditions to the ones prevailing at high redshift. 
As for the Press-Schechter theory in cosmology, the present theory provides a complementary 
approach and useful guidance to numerical simulations exploring star formation, while making 
testable predictions. 
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Introduction 



Since the seminal work of Salpeter (jSalpeter 19551 ) and its empirical derivation of the stellar initial mass 



function (IMF), defined as the number density of stars per log-mass interval, dN/dlog M, tremendous effort 
has been devoted to the characterization of the IMF in various environments, including the Galactic field, 
young clusters, star forming regions but also the Galactic bulge and halo or high-redshift galaxies, with 
the aim to identify the underlying physical mechanisms responsible for its behaviour (see e.g. Kroupa 2002, 
Chabrier 2003a, 2005 for recent reviews). A correct understanding of the IMF is one of the major, unresolved 
issues in astrophysics. Indeed, the IMF provides the essential link between stellar and galactic evolution and 
determines the chemical, light and baryonic content of the universe. On the other hand, various observations 
of the prestellar dense core mass function (CMF) have shown its striking similarity with the stellar IMF, the 
former one being shifted by a factor of ~ 2-4 towards higher masses with respect to the latter one (Motte et 
al. 1998, Testi & Sargent 1998, Johnstone et al. 2000, Andre et al. 2007, Alfves et al. 2007). This strongly 
suggests that the IMF might already be determined by the CMF, at the prestellar stage, and thus would 
be due to the dynamics of the parent molecular gas. In the present paper, we adopt this observationally 
supported perspective and we refer without distinction to the CMF/IMF for the derivation of the general 
theory, except in §7.1 where we specifically address this issue. 



1.1. Previous works 

Various attempts, by many authors, have been made to derive a general theory for the IMF, with the aim 
to understand its origin and thus to make robust predictions for the stellar mass distribution under various 
conditions. These theories are based on analytical or numerical studies and, without being exhaustive, invoke 
either gravitational fragmentation or accretion (Silk 1995, Inutsuka 2001, Basu & Jones 2004, Bate & Bonnell 
2005), turbulence (Padoan et al. 1997, Padoan & Nordlund 2002, Tilley & Pudritz 2004, Ballesteros-Paredes 
et al. 2006), purely independent stochastic processes (Larson 1973, Zinnecker 1984, Elmegreen 1997) or 
outflows (Adams & Fatuzzo 1995) as the dominant mechanism responsible for the CMF/IMF. 

The CMF/IMF has been inferred from hydrodynamical simulations using either the SPH technique 
and sink particles (Bate & Bonnell 2005) or grid codes and clump finding algorithms (Padoan & Nordlund 
2002, Padoan et al. 2007, Tilley & Pudritz 2004, Li et al. 2004, Ballesteros-Paredes et al. 2006), but with 
different conclusions. Padoan et al. (2007) conclude that the CMF/IMF produced in purely hydrodynamical 
simulations is too stiff, compared with the Salpeter value, whereas Tilley & Pudritz (2004), performing a 
Virial analysis, reproduce the Salpeter slope and Ballesteros-Paredes et al. (2006) find that the CMF/IMF 
depends on the Mach number. Padoan et al. (2007) also conclude that simulations including the presence of 
a weak magnetic field yield the proper CMF/IMF and thus that the presence of a magnetic field is mandatory 
to recover the correct Salpeter high-mass tail. 

Analytically, the most complete and successful approach to the derivation of the CMF/IMF has been 
proposed by Padoan & Nordlund (Padoan et al. 1997, Padoan & Nordlund 1999, 2002). It entails the 
following basic assumptions: protostellar cores result from overdensities due to shocks, and fluctuations 
of the underlying velocity field follow a Kolmogorov-like scale dependence. This enables them to obtain 
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the mass of the cores as a function of the scale. Then, these authors invoke the fact that the number of 
cores should depend on the scale L as M{L) oc L~^. Combining this latter relation with their mass scale- 
dependence, they obtain a distribution of the number of objects per mass interval, M = dN/dM, which 
follows a powerlaw. As a last step, they combine this distribution with the distribution of Jeans masses 
obtained from the probability density function (PDF) of the field. One of their main conclusions is that in 
the hydro dynamical case, the slope of the CMF/IMF is too stiff compared with the Salpeter value, while this 
latter is recovered in the weakly magnetized case, due to different magnetized shock jump relations. This 
approach, although retaining undoubtedly some interesting concepts, presents various shortcomings. First 
of all, the M{L) cx relation is not rigorously justified (see Elmegreen 2007). Furthermore, there is no 
proper integration over a fluctuation spectrum and the Jeans mass distribution should not be independent of 
the core mass inferred from the shock conditions. The shock conditions themselves are not well justified for 
the magnetic case since they assume that the magnetic field is always perpendicular to the shock. Finally, 
the result of the CMF/IMF being too stiff in the hydrodynamical case is in contradiction with the numerical 
simulations of Bate & Bonnell (2005) and Tilley & Pudritz (2004). At last, in the Padoan & Nordlund 
model, the turbulence is supposed to be taken into account implicitly in the shock velocity dependence and 
through the lognormal distribution of Jeans masses, but the turbulent support is not explicitly accounted 
for in this theory. 

In conclusion, none of the presently available analytical theories of the mass function provides a rigorous 
and robust foundation for the CMF/IMF. 

1.2. Aim of the paper 

Such a robust theory has been derived in the cosmological context by Press and Schechter (1974, 
hereafter PS), who have derived an analytical formulation of the number density of collapsing dark matter 
halos in the so-called hierarchical model of structure formation. In spite of its relative simplicity, this 
formalism has been shown to reproduce remarkably well the results of large computer simulations and has 
been applied successfully to various cosmological problems aimed at understanding galaxy formation. This 
formalism also provides precious guidance to observations or numerical simulations of structure formation 
in the universe. The first attempt to apply the PS formalism to the field of star formation was made by 
Inutsuka (2001). In the present paper, we extend this formalism to the characterization of the CMF/IMF 
of stcillar cores and of CO clumps, in the general picture of the gravo-thermal or gravo-turbulent collapse of 
a molecular cloud, with a selection criterion based on the thermal or turbulent Jeans mass for the stellar 
CMF/IMF, and on a simple density threshold for the CO clumps. Our anal}d;ical theory reproduces well 
the characteristic IMF derived from observations over the entire mass range from the high mass to the 
brown dwarf domains, with the proper characteristic mass and Salpeter high-mass power-law exponent. We 
demonstrate that the CMF/IMF is obtained from the statistical selection involving hydrodynamical thermal 
or non-thermal motions, characterized by a lognormal density distribution, and gravitational instability, 
through the threshold condition for gravitational collapse, which sets up the characteristic mass. This 
interplay between hydrodynamics and gravity leads to the transition from the power-law tail at large masses 
to a lognormal form at small masses, as found in the IMF inferred from observations (Chabrier 2003a, 2005). 
In the presence of large-scale turbulent motions, the slope of the power-law tail is determined by the spectral 
index of the turbulent power spectrum, and the proper Salpeter exponent corresponds to the present-day 
observed value in molecular clouds. The characteristic mass and the variance of the lognormal part are 
entirely determined by the local conditions prevailing in the cloud, namely the temperature, density and 
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thcrmal or turbulent rms velocity. Although based on the same general concept as the one suggested by 
Padoan & Nordlund (see references above), the present theory is more general and identifies shortcomings in 
these author approach. The paper is organized as follows: in §2, we derive the framework of our statistical 
formalism; in §3, we consider the simple density threshold and calculate the related mass spectrum of CO 
clumps; in §4, wc introduce the Jeans- type selection criterion. The general CMF/IMF analytical formulation 
is derived in §5, while the results are presented in §6. Section 7 is devoted to comparison with observations 
and with previous works. Section 8 concludes the paper. In order to facilitate the reading of the paper, the 
most important symbols and notations used in the paper are given in Table 1. 



X 


index of the IMF dN/dlogM oc M"^ (Salpeter index: x = 1.35) 


n 


3D power spectrum index of the velocity field (Kolmogorov case: n = 11/3) 


n' 


3D power spectrum index of the log-density (log p) field 


V 


exponent of the velocity dispersion vs size relation, rj = (n — 3)/2 


Li 


injection scale or system size 


R 


smoothing scale 


AS 


Jeans length 


R 




M 


cloud muss 


M'j 


Jeans mass 


M 


M/M'j 


M* 


Mass at which thermal and turbulent support are comparable 


P 


cloud mean mass density 


5 


logarithmic contrast density log{p/p) 


N{M) 


N{M)dM = dN number of stars of mass between M and M + dM 


Mo 




M 


Mach number at the cloud scale 


M, 


Mach number at the Jeans scale 


a{R) 


Width of the density distribution at scale R 




Width of the density distribution at scale Lj 


Cs 


sound speed 


Va 


Alfven speed 




root mean square velocity, T^-ms = Vo(i?/lpc)'' 


Vo 


root mean square velocity at i? = Ipc 



2. Statistical description 

2.1. PDF from supersonic turbulence 

In the PS theory of structure formation, the structures are identified with over-densities in a random 
field of density fluctuations 

n'5) = ^=L=expf-^V (1) 



where 5 = p/ p — 1, p being the local mean density, and a is the standard deviation of the distribution. 
However, whereas the density of the primordial universe is known to be very uniform, the density field in 
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molecular clouds is highly non-uniform, which considerably increases the degree of complexity to derive an 
analytical formalism. On large scales ( > pc) , the spectral line widths of molecular transitions observed in 
molecular clouds indicate highly supersonic motions characterized by Mach numbers M.—iy'^)^!'^ jCs > 1, 
where Cg = {kT / ^mHY^"^ ~ 0.2 (T/IOK)^/^ kms^^ denotes the thermal sound speed and ^ is 

the mean molecular weight. 

The exact nature of the turbulence in molecular clouds is still a matter of debate. Numerically, it has 
been found that non self-gravitating supersonic isothermal turbulence leads to density fluctuations which are 
Gaussian at all scales in terms of the logarithm of the density, 5 — log(p/p)0, yielding a so-called lognormal 
distribution in density (Vazquez-Semadeni 1994, Padoan et al. 1997, Passot & Vazquez-Semadeni 1998, 
Ostriker et al. 2001, Kritsuk et al. 2007). More specifically, it has been inferred from numerical simulations 
that 

where 5 = — (Tq/2, from eq.([2]), and 

= ln(l + bM^) (3) 

with b « 0.25. In this expression, the Mach number M can be either a hydrodynamical or an Alfvenic Mach 
number. 

Whether these equations accurately describe the dynamics of molecular clouds remains an open issue. 
In our approach, we will nevertheless adopt these expressions for reference. In any event, our formalism and 
our results can easily be applied to other distributions. 



2.2. Scale dependence of the density PDF 

In the cosmological literature, the dependence of a on the scale R at which the density field is smoothed 
has been established to be 

/"OO ^ 

a\R) = / S\k)W^{R)d'k, (4) 
Jo 

where Wk is some window function to be specified (see e.g. Press & Schechter 1974, Bower 1991, Padman- 
abhan 1993) and S{k) is the density power spectrum. This expression simply states how the variance of the 
density field is related to its power spectrum. In the context of molecular clouds, a similar approach, namely 
the delta-variance, has been proposed by Stutzki et al. (1998) and Bensch et al. (2001). 

Assuming that S"^ is proportional to k~" and using a window function sharply truncated in the fc-space, 
Wk{R) — 0{R^^ — k), where 0{z) = 1 if z > and otherwise, one gets cr^{R) oc i?" This expression is 
consistent provided n' < 3, since otherwise the integral diverges at small k (Padmanabhan 1993). 

Unfortunately, in the context of supersonic turbulence, the PDF stated by eq.([2]) does not entail any 
scale dependence, despite the fundamental need to specify the scale at which the quantity is measured. This 



'^Here we adopt a different notation from the cosmological case since 5 is not the difference between the variable log(p/p) 
and its mean. 
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is likely due to the fact that in turbulence the fluctuations are dominated by the large scales and thus, as 
long as the measurement is done at a scale which is small compared to the injection length Li, the deviation 
from eq.([21) remains small. In the case of weakly compressible turbulence, the density field undergoes small 
fluctuations. Therefore, in that case, it is possible to refer to the cosmological results (Padmanabhan 1993) 
and to use the expression given by eq. ([4]) , as pointed out by Hennebelle & Audit (2007) . When the turbulence 
is subsonic, the power spectrum index of the density field is known to be close to the Kolmogorov index 
obtained for the velocity riowev spectrum in incompressible turbulence, i.e. n — 11/3 (e.g. Bayley et al. 
1992, Kim & Ryu 2005) B In that case, however, the integral will diverge at small fc, as mentioned above, 
due to the fact that the larger scales contain most of the energy. Therefore, the integral must be truncated 
at the smallest values of fc, A:niin ~ 27r/Li, so that eq.(|4]) becomes: 

a^{K)^ 5'^{k)Wl{R)Sk^ 5^{k)A^edk^C\l-[^] , (5) 

where C is some constant. Note that Hennebelle & Audit (2007) have numerically checked that this expres- 
sion is reasonable. 

In case of supersonic turbulence, the power spectrum of log(p) has been calculated by Beresnyak et 
al. (2005) and Kritsuk (priv. comm.) in isothermal hydrodynamical and MHD simulations. The derived 
exponent, n', turns out to be close to the n = 11/3 value obtained in incompressible turbulence for the 
velocity field. Note that this is not the case for the exponent of the density field power spectrum which 
becomes much smaller when the Mach number increases (Beresnyak et al. 2005, Kim & Ryu 2005). Strictly 
speaking, in the case of supersonic turbulence, the use of eq.® is questionable since it is valid only for a 
periodic function, which is unlikely to be the case for the density field within a spatially finite molecular 
cloud. However, more generally, cr(i?) must obey the following properties. First of all, cr — s- when K — > L^; 
second of all, a must tend toward the expression given by eq.((5]) when ^ 0, implying C = CTq; third, 
CT (To when R ^ Li. 

Since a{R) is rather constrained and since the expression stated by eq.dH with C = CTq obviously satisfies 
all these constraints, it seems reasonable to assume that a is still given by eq.© even in the supersonic case. 
In any event, as will be shown later, the exact dependence of a upon R is not a crucial issue for the 
derivation of the CMF/IMF, as long as the aforementioned third condition is satisfied. Note also that our 
results regarding the mass spectrum of structures defined by a density threshold will remain unchanged as 
long as cr^(i?) = (Tq X /(I — {R/Li)^ ^^), with / any positive and monotonically increasing function such 
that /(O) = 0. 

Therefore, in the following, we consider the random field of density fiuctuations as given by: 

where a is given by eq. ([5]) with C = ctq. 

In our approach, the star-forming clumps issued from these large-scale turbulent motions are thus 



^Note that, in the present paper, the indexes n' and n of the power spectrum of the log-density and velocity fields, respectively, 
refer to the 3D value of the power spectrum. The usual ID Kolmogorov and Burgers values thus correspond to nio = n — 2 = 5/3 
and 4 — 2 = 2, respectively 
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identified with over-densities S = \og{p/p). The mass associated with the surdensity at scale R is 

M R^p^ R^pe\ (7) 

so that R = (M/p)^/^ exp(— 5/3). Note that since we use the sharp k-space filter, the relationship between 
R and M is ambiguous up to factor of unity. It is important to note that, unlike in the PS formalism, the 
mass depends not only on the scale R but also on the variable 5, a consequence of the lognormal, instead of 
uniform, underlying density distribution. 



3. Mass spectrum of structures defined by a density threshold 

Before we investigate the more complex case of self-gravitating objects, we first consider objects defined 
as in PS by a simple density threshold, p^. In the PS formalism, the structures above this threshold will 
have collapsed at the time of interest and will have formed a gravitationally bound object. In our case, 
this is different since in the ISM the gas is supported against gravitational collapse by various sources (see 
next section). However, structures above some density threshold may either undergo thermal instability, as 
studied by Hennebelle & Audit (2007), or simply be revealed observationally because they have a sufficient 
abundance of CO molecules. In any case, it is interesting to investigate this case because of its simplicity 
and because there are now numerous observations of the CO mass spectrum. 

The total mass. A/tot (-R), of the gas which at scale R has a density larger than pc = pexp{Sc) is simply: 

Mtot(i?) = / pexp(5) VRiS) dS, (8) 

where Li is the size of the system assumed to be comparable to the injection scale. The value Ssup is 
not consequential on the results and will be assumed to be equal to infinity, Ssup ^ oo, in the rest of the 
calculations. This total mass, Aftot(^), represents the mass of the structures of size larger than or equal to 
R and of mass larger than or equal to 

= Pci?'. (9) 

We are interested in counting the structures having a mass equal to il/]^. Some of these structures, 
however, have a non-zero probability to be included in larger structures which exceed the density criterion, 
and thus the number of these structures cannot be straightforwardly obtained by requiring that the total 
mass they contain be equal t o Mt^t- This is similar to the so-called cloud in cloud problem in the PS theory 
(|Bond 199lUjedamzik lOoih . We follow the approach of Jedamzik (1995) to handle this problem (see also 



Yano et al. 1996 and Nagashima 2001). 

Let M{M')dM' be the number-density of isolated structures of mass between M' and M' -t- rfM', which 
satisfy the aforementioned density criterion. More precisely, these structures are the ones which at scale 
R' ~ {M' / pcY^^ would have a density p ~ pc- The mass contained in such structures is M'Af{M')dM' . At 
scale R < R' , these structures contain regions whose density can be larger or smaller than pc- 

Let P{M'j^, M') be the mass fraction of a structure of mass M' which, at scale i?, has a density above 
the critical density. The mass of the gas having a density larger than pc at scale R and yet contained in a 
structure of mass M' as defined above is thus: M'J^{M')P{M'^, M')dM' . This yields a second expression 
for Mtot(i?): 

Mtot(i?) = i-/ M' N{M')P{M%,M')dM'. (10) 
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As emphasized by Jedamzik (1995), the exact value of P{M,M') depends on the choice of the window 
function and to some extent on the exact definition of a structure (see Yano et al. 1996). Jedamzik explored 
the influence of P{M, M') and found that its value has some limited influence on the result. However, in the 
case of a window function sharply truncated in the k-space, Yano et al. (1996) argue that P{M, M') — 1/2. 
This stems from the fact that, in that case, the smaller region inside the larger region has an equal probability 
to be overdense or underdense than 5c- Note that in the cosmological case, the Bond et al. (1991) result is 
exactly recovered if one assumes P{M,M') = 1/2. Since, as demonstrated in Appendix [C| it appears that 
the calculation of Yano et al. (1996) is also applicable in our case, we take P(Af|i., M') — 1/2. 

Equating eqs. ([8]) and (fTO)) and deriving the expression with respect to R, we obtain: 



mi'n) = -^4S^^% I AiS.mm. (12) 



This yields the following expression for the seeked mass distribution at scale R 

2p dR 1 da r 
M^dMf^y^a^dRJs^ 

where 

cr2 4 y y 20-2 

and, with the expression adopted in ^2.21 

da n'-ial ( R 



^(<5,i?)= (-1+ exp --3 + ---), (13) 



(14) 



dR 2a R \L,, 

The integral can be calculated analytically: 



A{S,R)dd= + — exp \ 2/ (^5) 



and is positive if 5c > —a'^/2 = 5, as in PS. As demonstrated in the appendix [Al when 5c < —a'^/2, which 
physically corresponds to a void instead of a cloud, 5c + a^/2 must be replaced by —{5c + a^/2) in eq. p^ . 
therefore insuring the positivity oi Af{M). 

With the aforementioned definition of (eqE]), we get for the mass spectrum of structures defined by 
the density threshold 5c: 

MiMi) = ^'"'-^"g (Mk)'^ X U+'l) oxp {-'''-f"\ . (16) 

(A/a^ 3v^a» \M„I \f,J V 2/ 1^ 2a' J ' ' ' 

where Mq = pLf is the whole mass contained within a volume and pc = pexp((5c). This expression is 
very similar to the one derived by PS in the context of a uniform density background. As in PS, we identify 
a power-law contribution and a Gaussian truncation around the threshold 5c, and we recover the fact that 
the fraction of bound objects of mass greater than Mr is proportional to 5c — 5. When the scale R ^ Li, 
then a and thus there is also a Gaussian cut-off for large-mass structures. Note that in the limit when 
5c ^ 5 = —(7^/2, JV{M) except if cr —> which precisely occurs when R — + Li. This implies that, when 
the threshold density 5c — > 5, all the mass lies within a unique structure of size Li, as expected. 
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On the other hand, there is no cut-off for smaU mass structures since they arise because of turbulent 
fluctuations which are scale free. This remains valid as long as the turbulent cascade remains self-similar, 
which implies that the scale R must be large with respect to the dissipation scale. Another important 
difference with the PS formalism is that time dependence is not taken into account. This point will be 
discussed in § 14.51 for the case of self-gravitating fluctuations. Here, we note that in the context of non 
self-gravitating structures defined by a simple density threshold, no time evolution is required to compute 
the mass spectrum of these structures. 

Equation naturally yields the scaling relation A/'( A/) cx Ad^'^ with (3 — —2 + {n' — 3)/3, so the power- 
law exponent is now affected by the spectral index characteristic of the logarithmic density power spectrum. 
This exponent is identical to the exponent obtained analytically by Hennebelle & Audit (2007, their eq.(15)) 
for subsonic turbulence and is in good agreement with the mass spectrum inferred from numerical simulations 
(Hennebelle & Audit 2007, Hennebelle et al. 2008, Heitsch et al. 2008). The exponent of the logarithmic 
density power spectrum has received only little attention but, as mentioned previously, values around the 
Kolmogorov exponent, n = 11/3, have been inferred from numerical simulations (Berezniak & Lazarian 
2005, Kritsuk, priv. comm.). For n' = 11/3, we get f3 = —16/9 = —1.777. This value is remarkably close to 
the slope of the mass spectrum inferred for the CO clumps (Blitz 1993, Heithausen et al. 1998, Kramer et 
al. 1998). This indicates that the CO clumps very likely have a turbulent origin. Note that in Hennebelle 
& Audit (2007), it was proposed that fluctuations of the warm neutral medium (WNM) induced by a sub- 
transonic turbulence are amplified by thermal instability once the density reaches the instability threshold. 
The present calculations indicate that the two mechanisms, namely subsonic turbulence followed by thermal 
instability, and supersonic turbulence lead to similar mass spectra. It is therefore difficult from observations 
of the mass spectrum only to determine which process is dominant. Note that the two mechanisms are 
actually not exclusive from each other. 



4. Selection criterion. Thermal and turbulent Jeans mass 

Unlike in the cosmological case, where the gas is very cold and not turbulent, the gas in the interstellar 
medium is supported by a combination of thermal pressure, turbulence and magnetic field. In this section 
we examine these various supports and we derive criteria to be fulfilled by the structures in order to collapse. 



4.1. Thermal Jeans mass 

In order for a cloud to collapse, its mass must be larger than the thermal Jeans mass 

where 



Mj = aj^^ = «J^= exp ( - 2 ) = ( " 2 ' ' ^^^^ 



MO ^ aj -3= « 1-0 aj {-^f^ (-^)-i/2 ( " Mo (18) 

where aj is a dimensionless parameter of order unity which takes into account the geometrical factor. In the 
absence of turbulent support, the clumps with mass larger than a Jeans mass will eventually collapse and 
form gravitationally bound objects. This implies a lower limit on the local density fluctuation, M > Mj — 

M° exp(-J/2) ^ ^ > -2 In {M/M°) . 
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Therefore, using eq.([7|), we obtain two equivalent conditions for the star forming collapsing structures: 
M>M^ = af^i? and ^ > ^ _2 In (^^) (19) 

where X°j = a]/^Cs/VGp « 0.1 (T/IO K)i/2 (^/2.33)-i/2 (n/10* cni-^)-!/^ pc is the thermal Jeans 
length 1. 

An important difference with the case of §3 is that the threshold now depends on the scale R. In 
particular — > oo as i? — *■ 0. Physically, this means that it is more difficult to have a gravitationally 
unstable object at small scales because of the thermal support. 



4.2. Turbulent Jeans mass 

If turbulence is significant, the turbulent support must be taken into account. As Tilley & Pudritz 
(2004), we use the Virial theorem to decide whether or not a structure is going to collapse. Neglecting the 
surface terms, the Virial theorem can be written as: 



2 dt 



2 _ 2Ecin + Epot + (3Pth + £;mag)V (20) 



where V is the volume and i?cin includes the contribution from turbulence, ^ 1/2 p{V^g), where (Kms)^^^ 
is the turbulent rms velocity, ii^pot is the gravitational energy and E'mag the magnetic energy. It is usually 
admitted that the structures which collapse are the ones having (Pl/dt^ < 0, although this criterion is not 
entirely rigorous, in particular if we were to apply it to one specific object. It is, however, a reasonable 
approximation for a statistical approach over a large population of objects, as in the present context (see 
e.g. Dib et al. 2007). This leads to 

(KL>+3(C.)2<-i?pot/M (21) 

Therefore, within our analytical formulation and our statistical description of the IMF, turbulent support 
can be included under the usual form of an effective sound speed (or equivalently effective pressure) (Chan- 
drasekhar 1951, Bonazzola et al. 1987, Vazquez-Semadeni & Gazol 1995 j^, 

C,,eff = [(C.)2 + i(y,L)]i/^ (22) 

The turbulent rms velocity is observed to follow a power-law correlation with the size of the region, the 
so-called Larson-type relations 



27J 



{VrD = Vo' ^[y^^) (23) 



with Vq ~ 1 kms ^ and rj ~ 0.4-0.5 ([Larson 19811 ). Note that, strictly speaking, the Larson relations are 



representative of the molecular gas at the whole cloud scale. Given the low efficiency of star formation and 



^ Strictly speaking the thermal Jeans length is ^/nCs / \fGp 

^Note that, because of the anisotropic nature of turbulent support, the effective sound speed formulation should represent 
an upper limit of the true non-thermal contribution. 
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the relatively high densities at which star formation is observed to take place, n ~ 10^-10^ gcm~'^ (Motte et 
al. 1998, Andre et al. 2007), it is not clear whether the Larson relations are representative of star formation 
regions. Indeed, the observed line-width vs size relation for prcstellar massive cores in star forming regions 
show slightly higher densities and velocity dispersions than predicted by the Larson relations (Caselli & 
Myers 1995). Therefore, the use of the Larson relations, although reasonable, should be considered as simply 
indicative. 

Using the fact that (^^^^^(i?)) = ^!^ij^Pv{k)dk cx where Pv{k) oc fc^" is the velocity power 

spectrum, and eq. , the exponent ry is related to the (3D) index of turbulence n by: 



(24) 



Note that, strickly speaking, the index of turbulence n which appears in this expression, which is related to 
the power spectrum of v, is not necessarily the same as the one previously introduced, n', which is related 
to the less standard power spectrum of log p. However, as mentioned earlier, numerical simulations seem to 
find that both indexes are rather similar. 

The aforementioned values of r/ thus lie between the ones corresponding to a Kolmogorov, 77 — 1/3, and 
a Burgers, 77 = 1/2, value, pointing to mildly to highly supersonic conditions in star forming clouds. Recent 
high resolution simulations of non-magnetized isothermal supersonic turbulence (Kritsuk et al. 2007) yield 
■q ^ 0.4-0.45 {n ~ 3.8-3.9). In the Kmit where thermal support can be neglected (valid for massive stars), we 
obtain a turbulent Jeans mass 

^3 / D \ S') 



tuib 



aj- 



^0 



(R_Y 



(25) 



^33G3pexp((5) Vlpc/ 

As for the pure thermal case, the structures which collapse in the turbulent case are the ones such that 
M > Mj^turb, which now yields: 



M > m: 



R 



2/3 

' 3G 



R 
Ipc 



2-n 



R 



and 



5>5'r. 



In 



2/3 



1/2 



( — ) 

3 Gpi?2 4pc^ 



2») 



(26) 



4.3. General case 



In the general case where both thermal and turbulent supports contribute, the condition for collapse, 
M > Mj, becomes 



M > aj 



{C,r + iVi/S)iR/lpc)'^ 



3/2 



v/G3pexp(^) 



which, with cq.([7]), implies 



2»7 



R 



(27) 



(28) 



and 



c^cc_i ( aTiiCsr + iVo'm (fl/ipc)^"] 



(29) 
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An important point to be stressed is that, in our approach, we select the regions of the gas at the very 
early stages of star formation, which will collapse in the future because the combination of all supports is 
not sufficient to balance gravity. In particular, the dense cores themselves, which are observed not to be 
very turbulent, represent in our approach the (collapsing) evolution from a selected initially more dilute and 
turbulent piece of fluid. This piece of fluid, which fulfills our criteria for collapse, is properly accounted for 
in the present theory, while the observed prestellar cores represent a more evolved state, of which physical 
properties are not described in the present theory. Similarly, it is known that turbulence dissipates in about 
one crossing time, so one may worry about not taking this effect into account in the theory. This can be 
understood as follows: if a piece of fluid contains too much turbulence, it will expand and gravity will be 
unable to take over. On the other hand, if the turbulence is not sufficient, gravitational collapse will proceed 
eventually. The decay of turbulence does not affect this selection process, unless perhaps for the cases where 
turbulence is just sufficient to balance gravity. 



4.4. Magnetic field 

The presence of a magnetic field, known to be dynamically significant in star formation (Crutcher 
1999), will modify to some extent this general picture of the origin of the CMF/IMF. Taking into account 
the magnetic field in such a theory requires to know exactly how magnetic field and density correlate. In 
numerical simulations, a broad correlation has been found between density and magnetic field, yielding 
B (X y/p (Passot et al. 1995, Padoan & Nordlund 1999). Observationally, using Crutcher's (1999) data, Basu 
(2000) has shown that the magnetic intensity follows very closely the relation B oc p^^^AV, where AV^ = 
(Kins)^^^ is the velocity dispersion. This suggests that the Alfvenic Mach number, Ma — AV/Va, where 
Va ~ B / [AirpY/"^ is the Alfven velocity, may be approximately constant in molecular clouds. Basu proposes 
an explanation based on magnetic flux and mass conservation and approximate mechanical equilibrium 
along the field lines. His calculation shows that the Alfven velocity can be written as Va = (YaI^s) x 
{CsY + Krns/3- Application of the Virial theorem (eq |20|) yields 

3 Cl + (KL) + ^ < -E^ot/M. (30) 

Therefore, in order to take into account the magnetic support in expression (j27p . we just have to replace 
[[CsY + Kms/3] by [{CsY + + Assuming the aforementioned dependence for the magnetic 

field, inclusion of magnetic support in our theory thus simply leads to an expression for the Jeans length and 
the Jeans mass which is proportional to the one derived in the hydrodynamical case. All the calculations 
conducted in this paper could thus, in principle, be generalized to the MHD case by, e.g., simply replacing 
aj by aj X [1 + (V2/Cs)^/6]'^/^ or, more generally, by simply rescaling the sound speed and the rms velocity 
by different factors. 

The magnetic field also modifies the value of the Mach number which appears in the width of the 
lognormal distribution. This implies that, in this simple approach, the magnetic field, because of the magnetic 
support, reduces the width of the density PDF (Vazquez- Semadeni et al. 2005). Although a more realistic 
treatment in which the magnetic field distribution would be properly taken into account is strongly needed, 
it seems difficult to perform it at this stage. Meanwhile, the present approach has the virtue of simplicity, 
while based on observational arguments. 
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4.5. Time dependence issue 

In the PS formalism, time dependence is accounted for by relating the collapse epoch of the perturbation 
to its density contrast. Although this somehow simplistic approach can be improved (see e.g. Audit et al. 
1997), it nevertheless yields satisfying results (Efstathiou et al. 1988, Lacey & Cole 1994). 

In principle, such a time dependence could also be obtained in our case by selecting at a time t not the 
gravitationally unstable pieces of fluid, but rather the ones that had time to collapse and form a singularity. 
However, the time evolution is more complex in molecular clouds than in the cosmological case, because 
of the various sources of support of the gas against gravitational collapse. Therefore, deriving the time at 
which collapse really occurs is a complicated task. Ignoring time dependence in our theory implies that the 
distribution of collapsing objects obtained at the very early stages by our threshold conditions is the one 
obtained once star formation is completed within the cloud of interest. 

A related problem concerns the dense core formation, subsequent accretion and possible merging (see 
Dib et al. 2007 and Peretto et al. 2007) not taken into account in a time independent analysis. As emphasized 
in § 14.31 the present analysis relies on a simple statistical counting of the smallest pieces of fluid which are 
dominated by gravity. Dense cores will form out of these pieces by progressively accreting the related initial 
reservoir of mass. Our analysis excludes any external accretion of gas which is not included in this initial 
reservoir. 

This lack of time dependence is certainly a limitation of the present theory. However, the agreement 
between this theory and the observed CMF/IMF, as shown below, seems to suggest that time evolution 
should not drastically affect the initial mass spectrum of collapsing structures. 

5. Mass spectrum of self-gravitating objects: derivation of the CMF/IMF 

5.1. Analytic formulation 

5.1.1. Physical assumptions 

Let us consider a region of scale R. The places where 5 is larger than the aforederived critical threshold 
(5|j, as defined by eqs. fTO)) . (pS)) and (P^)) . contain more than one Jeans mass and are going to form stars of 
mass smaller than or equal to Af^. 

This is because we assume that the final mass of the cores which form is equal to the mass associated 
with the clouds containing only one Jeans mass, since clouds which contain initially more than one Jeans 
mass are likely to fragment into several objects, whose number is more or less equal to the number of Jeans 
masses contained in the cloud. Therefore, all points which, at scale i?, have a density contrast larger than 
are going to form structures of mass sm.aller than or equal to Mj|. This can also be understood in the 
following way. Consider a cloud which at scale R contains about one Jeans mass, implying that its mean 
density at scale R is about MjjR? . At smaller scales, it may happen that the cloud is not uniform but is 
composed of smaller, denser cores embedded into a more diffuse envelope. If these denser cores contain one 
Jeans mass, the end product of the collapse is likely to be a cluster of objects whose mass is close to the 
mass of the smaller/denser cores and not to the mass of the object at scale R. We note here a fundamental 
difference with the cosmological case, where it is assumed that the mass of the final objects is equal to 
the mass of the biggest cloud which satisfies the appropriate conditions. This difference arises from the 
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fact that the density contrast between the standard molecular gas and the star itself is about 18-20 orders 
of magnitude. Note that in our approach, as mentioned earlier, we neglect any further accretion on the 
prestellar cores. More precisely, it is assumed that most of the mass which will eventually be accreted onto 
the prestellar dense core is contained within the initial Jeans mass reservoir. 

Given the aforederived threshold conditions for collapse, the places where S is smaller than (5^ contain 
less than one Jeans mass. Eventually, these regions can either form a mass bigger than A/^ or form no 
structure at all. Therefore, in the present context, what we are interested in is to find out the total mass, 
Mtot{R), which is going to form structures of mass or less. As just mentioned, this corresponds to the 
places with a density fluctuation S > 



The mass contained within structures of mass M < is equal to the mass of the gas which, smoothed 
at scale R, has a logarithmic density larger than S^. We have 



In a way similar to the approach followed in §3, the mass Aftot(^) can also be estimated by counting directly 
the self-gravitating clouds of mass smaller than Mf^. The number of such structures is Af{M')P{R, M')dM', 
where M{M')dM' is the density of structures of mass between M' and M' + dM' and P{R, M') is the 
probability to find a gravitationally unstable cloud of mass M' embedded inside a cloud of gas which at scale 
R has a logarithmic density larger than 5^. Therefore 



Note that the integration is from to M'j^ because if a self-gravitating cloud of mass M contains a self- 
gravitating cloud of mass M', then, as explained previously, we assume that an object (or few objects) of 
mass M' will form instead of an object of mass M . 

Estimating the probability M') is not straightforward. It seems possible to formulate it in a way 
similar to Jedamzik's (1995) formulation of P{M,M') in the cosmological case (his equations 8a and 8b). 
However, the complexity of the corresponding expression would lead to an equally complex result from which 
it would be difficult to extract the basic physical principles. On the other hand, Jedamzik finds significant 
deviation from the case P{M,M') = 1/2 only for structures 6 to 8 orders of magnitude smaller than M^, 
0. In the present case of interest, the IMF typically entails masses at most 3 orders of magnitude smaller 
than the mean Jeans mass. Therefore, it seems reasonable, in the present first step calculations, to follow 
the most simple approach. In the following, we thus assume that P{R,M') = 1. This means that we make 
the assumption that any self-gravitating cloud of mass smaller than is embedded into a cloud which, 
smoothed at scale R, is Jeans unstable. In other words, we assume that the Jeans unstable clouds are not 
isolated but are embedded into bigger Jeans unstable clouds (containing more than one Jeans mass). This 
assumption is further justified in Appendix [P] As will be shown below by our results, this seems to be a 
very reasonable assumption. 



5.1.2. Analytical expression 




(31) 




(32) 



^Interestingly enough, the discrepancy seems to decrease when his n decreases and in particular when it becomes negative, 
which is the case for turbulent fluctuations. 
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Equating eqs. |3T|) and ((32| and deriving the expression with respect to R, one gets: 



p dR 



(33) 



as the general expression for the number density distribution of star-forming coUapsing regions under the 
aforedefined selection criterion for collapse. 



As seen, two terms appear in eq. 



Since 6% is a decreasing function of R, the first term is positive. 



The second term is identical to the one which appears in eq. pTI) . As shown in appendix [Bj this term 
becomes significant only when R ~ Li, i.e. when the size of the structures becomes comparable to the 
size of the system itself. Clearly, in this limit the precise dependence of cr on i? becomes crucial and the 
statistical approach becomes rather questionable. From an even more fundamental point of view, this raises 
the question of the exact definition of the isolated fragmenting systems that are considered in this work, 
and how they are connected to the surrounding medium. We thus ignore this second term in the following, 
except in § 16.31 where it cannot be avoided. This implies that the spatial scale R must be small compared to 
Li or equivalenty that the mass of the structures must be small compared to the mass of the system itself. 



5.1.3. The normalization problem 

In order to check whether the expression we get for the mass spectrum of collapsing structures is correctly 
normalized, we consider a situation where the injection scale Li — > cxd. In this limit the system is infinite. 
Therefore, all its mass, pLf, will eventually collapse, because even a low density piece of fluid is contained 
within a Jeans mass. Then, we must have: 



POO 

pLi =Ll N{M)MdM , 
Jo 



which, with eq. ([31]) . leads to 



PLI = pL\ 



e,qy{5)VR{6)d5, 



(34) 



(35) 



It is easily shown that the integral is equal to 1 since it represents all the mass within the system. Thus, the 
equality is satisfied and our expression is adequately normalized. 



5.2. Mass spectrum with purely thermal support 

In case of pure thermal support, the mass spectrum of gravitationally bound objects derives from 
eqs. (HI]), (dH) and (|19p. and eq. (|55|) becomes, ignoring its second term: 



2-pM^j 



(M^)3 V2^CT 



exp 



[21n(Afj^/Mg)]^ 

2a2 8 



(36) 



which can also be rewritten 



- 16 - 



where M = M/M^. 

We clearly see from eq. ([37|) that the mass spectrum involves two contributions, namely a power-law 
with an index —3 and a lognormal contribution given by the term — 21n(Af)/cr^. The first contribution is 
dominant when <C M ^ M+, where 

M±=exp(±^a2), (38) 

whereas the latter one eventually becomes dominant both at very large {M M^) and at very small 
(M <C M~) masses, where it produces an exponential cut-off. The cut-off at small masses had been 
previously identified by Padoan et al. (1997). It arises from both the lognormal nature of the density 
distribution and the threshold condition for collapse (eqs fT9l [26l [29|) . as demonstrated in the present paper. 
This provides a rigorous foundation for this p eculiar form of the p ower-law exponent of the CMF/IMF, 



sometimes invoked empirically in the literature (jMiller fc Scalo 19791 ). 



This clearly demonstrates that the stellar CMF/IMF results from two contributions, a power-law which 
dominates in the aforementioned mass range, and a lognormal form, which becomes important at very small 
and very large masses, as characterized by the transition mass (eq i38|) . Equation ([37|l also shows that 
the mass spectrum of bound objects issued from a purely thermal collapse has a much steeper distribution at 
large masses, A/" cx M~^, than the one given by the Salpeter CMF/IMF, M{M) oc M~^-^^. It also highlights 
the importance of the characteristic scale of the system on the mass spectrum, through the scale-dependence 
of the variance a (see §2). Finally it also clearly shows the importance of the Mach number on the mass 
spectrum of collapsing prestellar cores. Indeed, with eqs.([3]) and ([38]) . we get: 

M± = (l + 5X2)±i, (39) 

Small-scale motions, i.e. small values of cr, will hardly produce any object far away from the mean Jeans 
mass. 



5.3. Mass spectrum with purely turbulent support 

In case the clumps are supported dominantly by turbulent motions, eqs. (fT^ . (fTS)) . pS)) and ([55)1 yield 



_ 2p (1 - v) ^6/(„-i) ^ exp(-aV8) 

(M0)2 (2?? + 1) * V2^a ' ^ ^ 

where M = M/Af^, M' = mI^'-^'^'^^M, ai = (1 + ri)/{2ri + 1), 03 = (ry - l)/(27/ + 1), as = 6/(2?? -I- 1), and 
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where the ratio (Vo/lpc)'' is given by the aforementioned Larson relation (eq |23l) . Roughly speaking, the 
effective Mach number measures the relative importance of turbulent versus thermal support contribu- 
tions at the Jeans length scale. Note that, according to the discussion of §4.4, this effective Mach number Al, 
can be renormalized to take into account the presence of a magnetic field. This requires, however, a proper 
knowledge of the exact dependence of both the uniform and fluctuating components of the magnetic field 
(or Alfven velocity) upon the thermal and non-thermal contributions of the velocity dispersion, respectively. 

From comparison between eas. ([57|) and we first note that the introduction of the turbulent con- 
tribution into the effective sound speed modifies the exponent of the power-law term, through the Larson 
exponent. Interestingly enough, the aforementioned favored value for supersonic turbulence, 77 — 0.4-0.45 
(Kritsuk et al. 2007) yields exactly the Salpeter coefficient, dN/dM cx M~^-^'^^\ Indeed, 1 + x = 3ai — 2.33- 
2.35, bracketed by the Burgers, 3ai — 2.25, and Kolmogorov, 3ai — 2.4, values. Expressed as a function of 
the (3D) index of turbulence, n, with the help of relation (|24p . the power law exponent of the mass spectrum 
as obtained in our calculations reads: 

(42) 



2n-4 

As for the thermal case, the precise value of the turnover mass, around which the CMF/IMF evolves 
from the power-law to the lognormal form, depends on the value of a and therefore on the Mach number 
M. Larger Mach values will produce larger numbers of small-scale collapsing clumps. 

The reason why the mass spectrum is stiffer when thermal support only is considered than when turbu- 
lent support is taken into account is simply because the turbulent support increases with the scale. Thus, a 
lot of intermediate to relatively large mass structures (of the order of or larger than the usual Jeans mass) 
which are unstable under purely thermal criteria are stabilized by turbulence. This support thus prevents 
fragmentation of these structures into several smaller structures, leading naturally to a shallower and broader 
mass spectrum in the high mass {M > Mj) domain. 



5.4. General case 



In the general case, where both thermal and non-thermal supports contribute, the mass spectrum now 
reads, from eqs. (fT2|) . (fTSl) . ([29]) . (f28|) and eq. (f33|) where, as explained previously, the second term has been 
dropped: 



Af{M) ~ 2 No 



1 



1 



1 + (1 - ri)Mlli^'^ 

I3 1 -I- (2?7 + 1)A12^2^ (l+A^2^2„)3/2 ^ 



X exp 



In 



exp(-CT2/8) 



V 



2a2 



27rcr 



where R = R/\^j, M = M/AfJ = i? (1 -f MlR^'^), = In <{ (1 + MlR^'^)/R^ \ and Mo = pl{M'}f 
This expression can also be rewritten: 

UiM) = 2No^ ^Hl-,)MlRy (M\ ^ - ^,exp(-.2/3) 

i?6 [1 + (27? + 1)M2^2„] I ^3 



27rCT 



(43) 



(44) 
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We see that the transition between the thermal and the turbulent dominated regimes, occurs when the 
radius R ~ M7^^^ and thus for masses around 

M* ~ 2{M,)-^^''. (45) 

For ~ 2, we get M* ~ 0.8 — 1. At masses larger than M*, we recover the power-law behaviour 
characteristic of the turbulent collapse, with the proper Salpeter value. This is easily verified for the case 
rj = 0.5, which yields Af{M)j^^^ cx M"^/^. The small-mass limit (i.e. M < M*), on the other hand, 
resembles the one of the purely thermal case. Therefore, at least for values of of the order of unity, we 
expect the mass spectrum of collapsing structures in the general case to be bracketed by the turbulent and 
thermal behaviours at large and small scales, respectively. 

6. Results 

In this section, we study eq. (|44p and its dependence upon the parameters, M (which enters the expres- 
sion of cr) and A^^. Recalling that the Mach number A4 is the ratio of velocity dispersion over sound speed 
at the scale of the whole cloud, whereas A^* is the ratio of velocity dispersion over sound speed at the scale 
of the Jeans length, we see that both parameters depend on the velocity dispersion. The first one increases 
with the size of the cloud while the second one increases with the size of the Jeans length. In order to 
investigate their respective influence on the CMF/IMF, we will first vary one of these two parameters while 
keeping the other one constant. Physically speaking, this corresponds to either considering clouds of fixed 
density but of various sizes {A4 varies but not M*), or changing the density but not the size (A^* varies but 
not Ai). For a sake of simplicity, we assume in these 2 cases that a — aQ. 

Finally, in order to investigate the influence of turbulence on star formation, we will vary simultaneously 
A4 and A^^ while keeping their ratio constant. This latter case corresponds to a cloud of fixed size and density 
but of varying velocity dispersion. In this case, the dependence of a{R) on R is properly taken into account. 

The results are expressed in terms of Mj as defined by eq. (fT8|) . Typically, for a cloud of density n ~ 10'^ 
cm~'^, Mj is about 3 Mq. For clouds of size ~ 1 — 2 pc following the Larson laws, typical Mach number 
values are Al ~ 6. Typical values of are around 1-2. These numbers must be considered as indicative, 
given the large uncertainties that remain on the exact conditions for star formation. In the following, the 
value of rj is taken to be equal to 0.4. 

6.1. Thermal and turbulent support 

Figure [T] portrays the mass spectrum obtained in the general case (eq. (|44l) '). for different values of Al*, 
namely = 0.1, 2 and 5, at a fixed Mach number M ~ 6. 

For the first value of Al*, the influence of the turbulent support is negligible (M* ~ 30) and the collapse 
is almost purely thermal. We see the steep cut-off at large masses due to the a = —3 exponent in the power- 
law, as well as the sharp exponential cut-off at small masses. As mentioned above, this latter term derives 
from the general formalism and from the critical density selection criterion, and thus is not specifically due 
to turbulence, but to the lognormal distribution characteristic of the density field. 

As expected, for Ail = 2 the behaviour at small masses is almost identical to the purely thermal case. 
The CMF/IMF peaks at almost the same position, although less intermediate mass stars form because of 
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Fig. 1. — Mass spectrum for = 6 and various values of A4^. For reference, power-law distributions 
dN/dM oc M~('^+^) are also plotted for a; + 1 = 2 and x + 1 = 2.3, about the Salpeter value. 
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the turbulent support, which prevents the collapse of large structures into smaller ones. At large masses, 
the CMF/IMF is less stiff than for the pure thermal case and has an index close to the Salpeter value, as 
mentioned earlier. In the case Ai^ = 5, the number of intermediate mass stars is even smaller, but the peak 
of the CMF/IMF remains weakly affected (as expected from eqn.([55])). 

6.2. Influence of the Mach number Ai 

Figure [2] shows the CMF/IMF for Ml = 2 and various values of the Mach number M. We stress that, 
since A^, is maintained constant, increasing the Mach number is not equivalent to increasing the value of 
Vq, but rather the size of the cloud. 

Clearly, supersonic turbulence strongly enhances the collapse of small-scale structures while the number 
of large-scale structures does not change significantly (for these values oi Ai). This stems from the fact that 
the transition mass, M~ , defined by eg. ([55]) . is modified by the variance cr, so that the exponential cut-off 
occurs at much smaller masses, for a given hydrodynamic Mach number A^, than for the thermal case. 

This shows that, if large-scale motions in molecular clouds are strongly supersonic, they will produce 
a large population of brown dwarfs. It is the fact that small-scale structures become subsonic (eq l23| . and 
thus of dominantly thermal nature, that they hardly form gravitationally bound objects far away from the 
mean Jeans mass, and thus that brown dwarfs, although issued from the same general formation mechanism, 
are not formed as efficiently as (low-mass) stars around the mean Jeans mass. In the high-mass regime, as 
mentioned previously, turbulence yields a shallower power-law tail than for thermal collapse. Mach number 
values ~ 6 yield the correct Salpeter index already for M >M'j, for a power spectrum index n ~ 3.8 
(77 ~ 0.4), as inferred from recent high resolution simulations of supersonic turbulence (Kritsuk et al. 2007). 
As expected, the general case indeed lies between the turbulence dominated collapse behaviour at large 
masses and the thermal behaviour at small scales. Not surprisingly, the characteristic mass of the spectrum 
is fixed by the characteristic mass for collapse, affected to some extent by the Mach number (which enters 
the variance), illustrating the role played by gravity in this collapse gravo-hydrodynamical picture. 

6.3. Support versus turbulent triggering 

A key, yet unsettled issue in the modern picture of star formation is to clearly understand the overall 
effect of turbulence on this process. On one hand, turbulence creates density enhancements (obviously related 
to A^, as given by eql39|) which tend to favor gravitational instability. On the other hand, the turbulent 
support tends to stabilize the gas (the turbulent support is due to A4* in our theory). Figures [1] and [2] 
show how the CMF/IMF changes with Al and AI*. However, if one increases the level of turbulence in 
a cloud of fixed size and mass, both parameters should change, while their ratio should remain constant. 
Figure [3] displays the mass spectrum obtained for three values of Ai, namely Ai = 3, 6 and 12, while keeping 
M/M^ = V3{L^/X°j)'' constant and equal to 6/^2 (i.e. X°j/L, ~ 0.1), leading to Ml = 0.5, 2 and 8, 
respectively. As expected, we obtain a combination of the behaviours observed in Fig. [1] and [2] the smaller 
the values of M and Al, the narrower the mass range of star formation around the characteristic mass Mj. 
Conversely, large Mach numbers promote the formation of a larger number of low- mass and high mass stars 
but produce less stars around Mj. 

Another important issue is to determine whether turbulence is globally promoting or quenching star 
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Fig. 3. — Mass spectrum for various values of A4 and A4* for a constant ratio = 
correspond to 3, 6 and 12 from the narrower to the broader distributions, respectively. 



4.24. Values of M 
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Fig. 4. — Ratio of the mass within stars over the initial mass as a function of the Mach number A4 for 
various values of (see text) (decreasing from top to bottom). 
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formation. Numerically (e.g. MacLow & Klessen 2004), it has been found that the global effect of turbulence 
on star formation is negative, although the turbulent support is much more effective if the driving is imposed 
at small scales. To verify this conclusion, we have computed the total mass included in the gravitationally 
bound prestellar cores, Mtot — J Af{M)MdM , for various values of A4 and various values of A4/A4^,, namely 
6/-\/2, 6/a/3, 6/\/4 and 6/\/5, which correspond to X^j/Li ~ 0.1, 0.18, 0.25 and 0.33, respectively. For these 
calculations, we have included the previously neglected second term in ea. psp and thus use the complete 
relation for M{M). As explained in §5.1.31 the integral MM{M)dM is constant. In reality, however, the 
system has a finite mass and size, so that the integral must be truncated when i? ~ or, equivalently, when 
M ~ pLii- Thus, the value of the integral depends on the injection scale. As discussed previously, integrating 
until R — Li is questionable, and we thus stop the integration at i? = 2Lj/3 = icut- We have verified that 
using different values of L^ut yields the same trends and qualitatively similar results. Quantitatively speaking, 
however, the results obviously depend on the value of R at which the integration is stopped. 

Figure [5] portrays the results of this global star formation efficiency as a function of the Mach number 
Ai. The top curve corresponds to the largest value of M/Ai^,, i.e. the smallest value of X^j/Li, whereas 
the bottom curve corresponds to the opposite. We see that, for a fixed value of M/M^, the higher the 
Mach number the smaller the star formation efficiency. In the same vein, for a given Mach number, the 
higher A4^ the less efficient the star formation. These behaviours are in agreement with and provide a 
theoretical foundation to what has been inferred from numerical simulations. We stress that, quantitatively, 
these results depend on the choice of Lcut- They demonstrate, however, that turbulence largely decreases 
the efficiency of star formation. 

7. Discussion 
7.1. Comparison with observations 

7.1.1. The shape of the CMF/IMF 

Figure [5] compares our analytical CMF/IMF fea H^ for 3 values of the Mach number, namely Ai — 6, 
12 and 25, with the IMF representative of the Galactic field and young clusters (Chabrier 2003a). This latter 
reflects the so-called system IMF since, in general, present limitations on angular resolution do not allow to 
resolve these systems into individual objects. On the other hand, our calculations are representative of the 
early stages of the collapse and do not consider possible subsequent subfragmentation of the prestellar cores 
into multiple objects. The same is true for most of present numerical simulations of star formation. In order 
to facilitate the comparison between the theoretical and the observationally derived distributions, this latter 
has been shifted in order to match the position of the peak of the analytical distribution. We will come back to 
this point below. The general agreement between the two types of distributions is striking. In particular, the 
transition between the power-law tail and the lognormal form is very well reproduced. This clearly assesses 
the validity of such a composite functional form for the stellar IMF, whose physical foundation has been 
demonstrated in the present paper. At the low-mass end, the theoretical and observationally-derived mass 
spectra start to differ noticeably below some mass value which depends on the Mach number, as analyzed in 
the previous section. A value A4 = 6 leads to a deficit of both very massive stars and brown dwarfs compared 
with the observed distribution, because of the too restricted mass range between the high-mass and the low- 
mass exponential cut-offs. A value M = 12 improves the situation and yields a fairly reasonable agreement 
with the observational distribution, although still underestimating the number of very small ( < 10^^ M©) 
brown dwarfs. A value A1 = 25 yields a good agreement with the observationally derived mass spectrum 
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Fig. 5. — Comparison between the theoretical IMF/CMF, dN/dlogM, (solid line) obtained with M = 6 
(lower panel), A4 = 12 (middle panel), = 25 (upper panel) and A4l = 2 and the stellar/brown dwarf 
system initial mass function (dotted line) of Chabrier (2003a). The peak of the this latter IMF has been 
adjusted arbitrarily to the one of each theoretical mass function. 
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over basically all the observationally probed mass domain. Although values ~ 25 are substantially larger 
than the ones usually associated with typical star forming regions, A4 ^ 6-12, two important points need to 
be considered in these comparisons. 

First of all, the statistics of brown dwarf detections below ^ 10^^ Mq are still small and determi- 
nations of the brown dwarf densities remain hampered by both observational (magnitude-limited samples, 
spectral-type effective temperature conversions, etc..) and theoretical (mass-age-L-Te determinations) 
uncertainties. The observationally-derived IMF in this domain thus retains a large degree of uncertainty. 

Second of all, the relation between the variance of the PDF obtained by turbulence, which determines 
the width of the CMF/IMF, and the Mach number is given by eq.(l3|), with b = 0.25. This relationship 
is extracted from 3D simulations of isothermal turbulence (Padoan & Nordlund 2002). Simulations taking 
into account a detailed treatment of the thermal properties of the gas in the molecular cloud, however, lead 
to an equation of state softer than isothermal, P cx pO-^-O-S^ least up to densities of a few lO'^ cm~^ 
(Glover & MacLow 2007). We thus expect a larger dispersion of the PDF and thus a larger extension of 
the CMF towards small masses with a more compressible non-isothermal gas, for the same Mach number. 
Furthermore, recent numerical simulations of supersonic turbulence with compressible driven modes predict 
a significantly broader PDF, for the same rms Mach number, than simulations with rotational driven modes, 
like the ones of Padoan & Nordlund (W. Schmidt, priv. com.). These results suggest that the relation Q 
underestimates somehow the aforementioned value of b and thus the variance a for a given A4. Finally, 
self-gravity could also possibly broaden the density PDF, as supported by numerical simulations (see e.g. 
Vazquez-Semadeni et al. 2008), although the densest part of the density distribution in that case departs 
from a pure lognormal form. 

Consideration of these issues should improve the quantitative agreement between the present theory 
and observations. 



7.1.2. From the CMF to the IMF: the core-to-star formation efficiency. 

As mentioned above, the theoretical and observationally-derived mass functions, dM/dlogM, have 
been arbitrarily adjusted peak-to-peak. After the shape of the CMF/IMF, we now examine the agreement 
between the real positions of the peaks of the CMFs, which define the most probable stellar mass scale. 
The comparison with the observationally-derived IMF on a solar-mass scale, however, depends on the exact 
value of the Jeans mass which, in our approach, defines the characteristic mass for collapse. Given the fact 
that the observationally derived system IMF peaks at roughly 0.3Mq, this would imply that the Jeans mass 
should be about IM© for Al = 6 and - 3 Mq for M =25. 

On the other hand, observations of prestellar condensations, as identified in dust continuum surveys, 
show that the CMF peaks at a mass about three times larger than the observational IMF, even though this 
peak determination still suffers from large uncertainties (Motte et al. 1998, Alves et al. 2007, Andre et al. 
2007, 2008). Given the similarity of the shapes of the CMF and the IMF, this implies that the transformation 
of gravitationally bound prestellar condensations into genuine stellar or brown dwarf systems involves some 
uniform star formation efficiency factor, e = M-^/Mcore, observed to be in the range e ~ 30-50%. Clearly, 
the process (es) responsible for the evolution of the CMF into the IMF remains ill-determined, although 
magnetically-driven outflows, expected to produce a mass-independent star formation efficiency factor in the 
appropriate range (Matzner & McKee 2000), offer an appealing explanation. This implies that the above 
estimate for the Jeans mass value adequate to match the observed IMF must be further increased by a factor 
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of the order of - 3, yielding ~ 3 Mq for = 6 and Af^ ~ 10 M© for M ^ 25 



7.1.3. Physical conditions 

Equation ([T5)) shows that, for star forming core conditions /i = 2.33, T = 10 K, n = 10'^ cm~'^, the Jeans 
mass is about 3 Mq, whereas it is about IMq for n — 10^ cm^'^. Fluctuations around these temperature 
and density values, as well as geometrical factors, however, will affect this determination to some degree, 
moving around the peak of the theoretical mass spectrum on a solar-mass scale. This would be consistent 
with the IMF being determined at gas densities in this density range. However, as discussed in §4.4, the 
presence of a magnetic field increases the critical mass for collapse, moving the peak of our CMF, normalized 
to the standard thermal Jeans mass feq fTS)) . to larger masses. Estimates of the critical magnetized mass 
for collapse, as done in various seminal works and textbooks (e.g. Spitzer 1978, Lequeux 2005), predict 
values for the relevant magnetic Jeans mass, Af$ c>c B^/n^, of the order of ~ 10 Mq for n = 10^ cm^'^ and 
-6 = 10 fiG. This implies that the typical gas density to produce a magnetic Jeans mass compatible with a 
peak of the CMF around 1 Mq, as found observationally, should be higher than 10^ cm~'^, with a typical 
value of the magnetic field slightly larger than 10 /xG. This is consistent with the densities characteristic of 
star forming clumps, n ~ 10^— 10^ cm~^ (Andre et al. 2007, Motte et al. 2007). On the other hand, the nice 
correlation observed by Basu (2000) suggests that the magnetic field also depends on the velocity dispersion. 
According to the analysis of §4.4, this implies that only the uniform magnetic component should enter our 
Jeans mass determination. Unfortunately, the value of this component is not determined observationally. 
Clearly, more knowledge on the structure of the magnetic field and its correlation with density and velocity 
is required at this stage in order to determine more precisely the value of the magnetic Jeans mass and to 
make a precise comparison between the characteristic mass scales, thus peak positions, of the theoretical and 
observationally derived IMFs. 



7.1.4- Expected dependence of the peak position in the cloud mass 

One surprising observational evidence, is that the peak of the IMF is fairly constant from one region to 
another (e.g. Chabrier 2003a, 2005). 

Here, we investigate how the peak of the IMF depends on the mass of the cloud. Since, as emphasized 
earlier, the peak occurs in the thermally dominated regime, we simply take the derivative of A/" as given by 
ea.(|37)). with respect to M, which yields: 

^-- = ^^p(-r')^(TTiW^- ^''^ 

This indicates that, for clouds characterized by high Mach numbers, the peak of the core mass function is 
roughly proportional to A4^^^^. Combined with eq. ((23|) . this gives Afpcak oc Since, typically, the 

CO clump density varies with their size as p (x L^"', with a = 0.7— 1 (Larson 1981, Heithausen et al. 1998), 
we get 

Mpeak = M°Mpeak OC L(^3/2)r,+a/2 (47) 

For 7] — 0.5, a = 1 or for rj — 0.4, a = 0.7, we find the same value (— 3/2)?7 + a/2 — —0.25 while for r/ = 0.5, 
a = 0.7, (-3/2)77+0/2 = -0.4. Since typically M cx i^-a.s ^j^ggg objects, we find that A^pcak oc A/-"-i""-2. 
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Thus, changing the mass of the clump by a factor 10^ changes the location of the peak of the IMF by a 
factor less than ~ 2. This partial compensation between the comparable increasing and decreasing scale- 
dependence of the Mach number and the Jeans mass, respectively, may be one of the reasons why the peak 
of the IMF appears to be rather constant over a wide range of stellar cluster conditions (see also Elmegreen 
et al. 2008). 

7.2. Comparison with previous works 

Although based on the same underlying picture of hydrodynamical- or turbulence-driven fragmentation, 
the theory presented in this work is more general than the one proposed by Padoan & Nordlund (2002). 
First of all, the present theory does not invoke shock conditions whose validity might be questionable, since 
they assume that the magnetic field is perpendicular to the shock. Note that, under this assumption, the 
shock jump conditions imply that B oc p, which is not compatible with the observations. Second of all, 
the present theory does not assume a one-to-one correspondence between the probability distribution of the 
turbulent gas density and that of local Jeans masses, identified as collapsing protostars. Regions of the gas 
which will collapse into a gravitationally bound object are properly accounted for, as they must fulfill the 
Jeans mass criterion through the proper selection condition for the density threshold. Regions where density 
fluctuations fail to fulfill this criterion will not be included into the mass spectrum distribution of collapsing 
objects. Therefore, the theory insures a correct counting of the collapsing structures. Finally, turbulent or 
magneto-turbulent support enters explicitly in our theory whereas it is interesting to note that this is not 
the case in the Padoan & Nordlund theory. Note that their original CMF/IMF (Padoan et al. 1997) closely 
resembles our pure thermal case. It is very likely the reason why, in the pure hydro case, their high mass tail 
is dominated by the steep dN/dM oc M^^ power-law, yielding an underestimate of the number of massive 
stars relative to low-mass ones. 

As in the Padoan & Nordlund theory, we recover the fact that the CMF/IMF is shaped by the product 
of two contributions, namely a power-law form, which dominates at intermediate scales, and two exponential 
cut-offs, one below about the characteristic mass for collapse and one at large scales. However, we note 
that this low-mass cut-off is not the result of turbulence, but is the result of the Gaussian, more precisely 
lognormal distribution of the density field, and of the threshold condition for collapsing structures. Although 
turbulence does produce such a type of distribution, other mechanisms, like wave propagation or gravity, 
could presumably lead to similar distributions. Turbulence, however, leads to a significantly smaller value 
for the low-mass cut-off, compared with the thermal case, promoting the formation of low-mass objects, in 
particular brown dwarfs. We also note that in the Padoan & Nordlund theory, the high mass part of the IMF 
is a pure power law whose index is determined by the shock conditions, with x = 3/(6 — n)[f|, whereas our IMF 
contains two contributions. A pure power law, of index x — {n + l)/(2n — 4), and a lognormal contribution 
resulting in an exponential decrease at the injection scale, as stated by eqs. (|36l) and ((43|) . Interestingly, in 
the Padoan & Nordlund theory, x increases with n whereas in our case x decreases with n. This fundamental 
difference is due to the fact that in the Padoan & Nordlund theory, turbulence is always promoting star 
formation since turbulent support is not included, whereas in our case, turbulence is also supporting the 
gas against gravitational collapse, leading globally to a negative effect on star formation, as shown in §6.3. 
Testing these qualitatively different predictions on the n-dependence would bring precious information on 



®We recall that n is here the 3D value of the power spectrum index. It thus corresponds to n = /3 -f 2 in the Padoan & 
Nordlund theory, where x = 3/(4 /3). 
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the very role of turbulence in star formation. 

Our theory also explains why hydrodynamical simulations of turbulence do recover the correct Salpeter 
power-law slope at high masses (Tilley & Pudritz 2004, Bate & Bonnell 2005). It is also in agreement with 
the finding that the CMF/IMF depends on the Mach number (Ballesteros-Paredes et al. 2006). 

Finally, our approach provides a consistent explanation, within the same general framework, for the 
mass distribution of both self-gravitating structures and CO clumps. This fits well the recent results of Dib 
et al. (2008) where it is shown that clouds defined by a low density threshold present a flat mass spectrum 
of index ~ —1-7, whereas clouds having a higher degree of gravitational binding and defined by a higher 
density threshold have a stiffer mass spectrum. 

7.3. Restriction of the present work 

It is necessary to emphasize various limitations of our present theory. These issues require further 
detailed investigations. First of all, as known in the cosmological context, the window function used to 
perform the smoothing at scale R has some influence on the results (e.g. Nagashima 2001). The same is true 
for the probability P{M, M') of finding a structure of mass M inside a structure of mass M', as demonstrated 
by Jedamzik (1995), although Jedamzik finds a relatively modest influence of P{M, M') (although see Yano 
et al. 1996 and Nagashima 2001). Second of all, the turbulent Jeans mass is treated using a mean scale 
dependence. Ideally, this could be improved by considering a velocity distribution correlated to the density. 
The same is true for the magnetic field distribution. 

At last, another important aspect already mentioned in ^ 34. 5) ignored in this work as in most of the 
other theories of the IMF, is the time dependence issue in the star formation process. The present theory 
consists simply in counting the fluctuations of a given distribution, produced by an underlying density field. 
In principle, the density fluctuations should evolve with time and rejuvenate. In particular, one may wonder 
whether the small scales should not rejuvenate more rapidly than the larger one. We think that part of the 
answer comes from the assumption of ergoditicity, which states that spatial averaging should give the same 
results as time averaging. In particular, the ratio of small to large scale fluctuations should remain the same. 
This is, however, obviously not the case for the scales comparable to the size of the cloud itself. For those, 
a time-dependent theory seems unavoidable. 

In a related way, our approach does not consider any accretion from external sources nor any further 
fragmentation during the collapse. Both could make the IMF different from the CMF, although the magnetic 
field seems to reduce the fragmentation significantly (Hcnnebelle & Teyssier 2008, Machida et al. 2008, Price 
& Bate 2008). On the other hand, the observed strong similarity between the CMF and the system IMF 
suggests that this latter is determined essentially by the CMF and that its shape should not be drastically 
modified by accretion, other than for the matter already included in the core reservoir. As for subfragmen- 
tation of the cores into individual objects, it has been shown that taking into account the mass-dependent 
multiplicity frequencies observed in the solar neighborhood seems to provide the correct link between the 
IMF of unresolved systems and the one of resolved individual objects (Chabrier 2003a, 2003b). We expect the 
same to be true for our theoretical CMF/IMF, although identifying the exact physical processes responsible 
for this subfragmentation remains an open issue. 
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8. Conclusion 

In the present paper, we have derived an analytical theory for the stellar initial mass function and for the 
mass function of the CO clumps, based on an extension of the statistical Press-Schechter formalism derived 
in cosmology. Our theory provides a predictive theoretical foundation to understand the origin of the stellar 
IMF, and to infer its behaviour in various environments. The theory predicts that the CMF/IMF involves 
two contributions, namely a power-law tail and an exponential cut-off below about the mean thermal or 
turbulent Jeans mass, even in the absence of turbulence. Although thermal collapse produces too steep a 
slope compared with the Salpeter value, this latter is recovered exactly in the case of supersonic turbulence, 
for the appropriate observed or numerically determined power spectrum index. This corroborates the general 
gravo-turbulent picture of star formation, as initiated by Larson (1981) and developed more recently by 
Padoan & Nordlund (2002), where large-scale protostcUar clumps which contain several Jeans masses are 
dominated by supersonic turbulent motions and will fragment into prestcUar cores that produce the final 
stellar spectrum. The smaller clumps will have subsonic internal velocities, in agreement with the observed 
Larson's relations in molecular clouds, i.e. will be supported by thermal motions, leading to a turnover 
of the CMF/IMF about the characteristic effective Jeans mass. Turbulence favors the formation of both 
low-mass and high-mass structures but as a whole it has a negative effect on star formation, decreasing the 
overall star formation efficiency. We also suggest that the opposite, comparable scale-dependences of the 
Mach number and Jeans mass lead to a weak dependence of the location of the peak of the IMF upon clump 
masses, providing an appealing explanation for the observed rather universal behaviour of the IMF over a 
wide range of stellar cluster conditions. 

The success of the present theory in reproducing the CMF/IMF inferred from the observations of 
prestellar cores strongly suggests that the IMF is determined by the conditions prevailing in the cloud, 
temperature, density and scale-dependence of the velocity dispersion (i.e. charac;teristic velocity power 
spectrum), and thus is already imprinted at early stages in the cloud. The universality of the IMF, at least 
under present-day Galactic conditions, very likely arises from the universality of the self-similar nature of 
turbulence and from comparable characteristic cloud conditions, determined by the same dominant cooling 
processes. Effects like accretion, ejection, collisions, winds may affect to some extent the exact shape of the 
CMF/IMF, explaining possible statistical variations, but are unlikely to be dominant mechanisms. 
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A. Mass spectrum of voids 

Here, we show that the voids have the same mass spectrum as the structures. The voids are regions 
of gas where the density is smaller than the average density. The mass contained within structures of mass 
smaller than A//? = PcR^ is given by 

Mtot(-R) = / ' pexp(5) Vr{5) d6. (Al) 

J —CO 

while eq. (fTO|) remains applicable. Thus 



2p dR 1 da 

J\ \1V1R) = - 

From eg. p^ . we see that 



This implies that 



/oo 
A{5, R)d5 = 0. (A3) 
-oo 



<5c 

A{6, R)d5 = - I A{6, R)d6. (A4) 



Therefore, the mass spectrum of voids is given by 



B. The dVu/dR term 



Here, we investigate the role played by the second term which appears in ea. p3|) . From eqs. p2|) . p5)) . ([29]) 
and eg. ([55]) . eq. fiSl) should be: 
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where Li — Li/X'j and where, for simplicity, we have assumed that n = n' . 

We thus see that the second term becomes comparable to the first one only when {R/ Li)"~^xaQ/a^{R) ~ 
1, which happens when R~ Li, i.e. for structures whose size is comparable to the injection scale, expected 
to be the size of the whole system. 
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C. On the P(M, M') = 1/2 result 

Here, for completness, we essentially reproduce the Yano et al. (1996) result which shows that, for a 
window function sharply truncated in the k-space, one has P{M, M') = 1/2. 

Let us consider a subregion of the flow which, smoothed at scale R', has a density equal to pc- Let 
M' ~ pc{R'Y ■ The probability that a region of size R < R' has a density larger than and thus contains 
a mass larger than MJj, is given by: 

P(M]„ M') = exp I ( ^—^ ] \ d6, (CI) 



27r CTsub 




where asuh = (^^{R') - (J^(i?). 

Indeed, this shows that P{M^,M') = 1/2. 



D. On the P{R, M) = 1 assumption 

Here we give justifications of the P{R, M) = 1 assumption done in the paper. Let us consider a spherical 
cloud with a density profile p oc r^". The cloud mass thus grows with r as M{r) oc r'^^". 



The thermal Jeans mass, Mj, is proportional to 1/y^, and thus Mj oc r"/^. This implies that 

oc r(3/^)«-^ (Dl) 



M(r) 

Therefore, if a < 2, the Jeans mass grows with r less rapidly than M(r) and a Jeans unstable mass is thus 
embedded into a larger more gravitationally unstable cloud. This is precisely the assumption P{R, M) = 1. 

In case of turbulent support, Mj(r) oc (Cs^eff)^/^/P> ^'^d Mj(r) oc r^^+°'/'^. Thus, the Jeans mass grows 
less rapidly than M(r) if 3r/+(3/2)a-3 < 0, leading to a < 2-2??. With r] ~ 0.4-0.5, this yields a < 1-1.2. 

Since it seems reasonable to assume that most cloud density profiles should not be much stiffer than 
l/r^^~'^\ we conclude that our assumption P{R,M) = 1 is realistic. 
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